Evidence for hadronic deconfinement in p-p collisions at 1.8 TeV 
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Abstract 

We have measured deconfined hadronic volumes, 4.4 < V < 13.0 fm^, 
produced by a one dimensional (ID) expansion. These volumes are directly 
proportional to the charged particle pseudorapidity densities 6.75 < dNc/dr] < 
20.2. The hadronization temperature is T = 179.5 it 5 (syst) MeV. Using 
Bjorken's ID model, the hadronization energy density is ep = 1-10 it 0.26 
(stat) GeV/fm^ corresponding to an excitation of 24.8 it 6.2 (stat) quark- 
gluon degrees of freedom. 
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The observation of high total multiplicity, high transverse energy, non jet, isotropic 
events led Van Hove and Bjorken ^ to conclude that high energy density events are 
produced in high energy p-p collisions Q]. These events have a far greater cross section 
than jet production. In these events the transverse energy is proportional to the number of 
low transverse momentum particles. This basic correspondence can be explored over a wide 
range of the charged particle pseudorapidity density dN^/ drj in p-p collisions at center of mass 
energy ^/s = 1.8 TeV. The various measurements from the Fermilab quark-gluon plasma 
search experiment E-735 have already been published. In this letter, we present for the first 
time a coherent picture based on the relationship between the volume V, temperature T, 
energy density ep, and pions per fm^ ra^r emitted from V. Spectra of identified particles, vr, 
K, ip, p, p, A° A°, S~, E~ are used to extract the V, ep, and np values and to determine 
the strange quark content and relative yields of the hadrons. 

Previously the various individual measurements did not provide an overall understand- 
ing of these p-p collisions. Prompted by the new analysis of the initial collision, we have 
developed a self-consistent picture of hadronic deconfinement. This letter discusses: (1) 
The role of parton-parton (gluon) scattering; (2) The volume at decoupling, resulting from 
the one dimensional longitudinal expansion; (3) The number of pions per fm^ emitted by 
the source; (4) The hadronization temperature of the source; (5) The hadronization energy 
density of the source; (6) The number of quark-gluon degrees of freedom in the source; 
(7) The deconfined volumes and plasma lifetimes, estimates of initial energy densities and 
temperatures. 

Experiment E-735 |]5| was located at the C0 interaction region of the Fermi National 
Accelerator Laboratory (FNAL). The p-p interaction region was surrounded by a cylindrical 
drift chamber which in turn was covered by a single layer hodoscope including endcaps. This 
system measured the total charged particle multiplicity 10 < Nc < 200 in the pseudorapidity 
range |?7| < 3.25. A magnetic spectrometer with tracking chambers and time of flight 
counters, provided particle identified momenta spectra in the range 0.1 < pt < 1.5 GeV/c. 
The spectrometer covered —0.37 < 1] < -1-1.00 with Aip ~ 20° { (f is the azimuthal angle 
around the beam direction ). 

(1) Recently the E-735 collaboration has analyzed the charged particle multiplicity dis- 
tributions arising from p-p and p-p collisions over a range of center of mass energies 0.06 
< a/s < 1.8 TeV 0. Results at 1.8 TeV support the presence of double ((72) and triple (as) 
parton interactions. These processes increase the non-single diffraction cross section (NSD) 
from ~ 32 mb at ^/s = 0.06 TeV to ~ 48 mb at ^/s = 1.8 TeV. The variation of the double 
encounter and triple encounter cross sections (72 and with center of mass energy y/s is 
shown in Fig. 1. 
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FIG. 1. Comparison of the cross sections for single, double, and triple encounter collisions 
which increase cjnsd above 32 mb function of \/s. 

The multiplicity distribution is made up of three contributions corresponding to single, 
double and triple parton-parton collisions. Our work on multiparton interactions shows that 
the increase in the p-p inelastic cross section with energy is nearly completely accounted for 
by the increase in multiparton interactions. Previously this increase in the p-p inelastic 
cross section was ascribed to copious minijet production As the energy is increased, a 
decreasing fraction of the center of mass energy appears in the NSD part of the inelastic 
cross section. This may be due to the decrease of the Feynman x of the partons involved in 
these collisions. It is thus likely that gluons become more involved with increasing energy 
leading to rapid thermalization P]. 

(2) To measure the hadronization volume V , pion HBT (Hanbury Brown, Twiss) corre- 
lation measurements were made as a function of both P^tt = P1+P2 the total momentum of 
the pion pair and of dNd di]. The Ptttt momentum dependent results are shown in Table I . 
Rg is the Gaussian radius parallel to the beam, r the Gaussian lifetime, and A the chaoticity 
parameter. The lifetime r can be viewed as a measure of the radius perpendicular to the 
beam [^. The increase of Rq and r with decreasing P^^^ is the characteristic signature for 
the expansion of the pion source [|l^. The dependence of Rg and r on dNc/drj is shown in 
Table II 0. A clear increase of Rg with dNc/drj is evident. The dependence of Rg and r 
on Ptt^ and dNc/drj is consistent with a one dimensional (ID) longitudinal expansion of the 
pion source. The effect of a ID expansion on the Rose Einstein correlation has been calcu- 



lated for a massless relativistic ideal gas |T2[. This calculation provides correction factors 
in and ir to our values of Rg and r obtained from the HBT analysis. Both and ir are 
a function of P,r7r and Arj, where At] is the spectrometer aperture. The cylindrical volume 
V of the pion source is V = 'K{irTY 2 IrRg where Rg varies with dNc/drj and irT reaches 
an asymptotic value for the larger dN^/ drj values. From our data Rg = e + h dNc/ dt] where 
e = (0.0788 ±0.013) fm and h = (0.0730 ± 0.011) fm and x^/NDF = 3.09/4.00 as shown in 



Fig. 2. We neglect e since h dNc/dt] is 6 to 20 times larger than e. The cylindrical volume 
becomes, 



V 



IT 



ilr^ 2 iRh dNjdri 



(1) 



2.5 



7.5 - 



0.5 - 












20 



25 



30 



5 10 15 
dN^dri 

FIG. 2. Dependence of the Gaussian radius Rq on dNc/drj. The gluon diagram indicates that 
two gluons are required to form two pions. 



The largest measured value of r = 0.95 fm is used to evaluate V . We estimate the £ 
factors using the extrapolation procedure (Ptttt— *-0) outlined in Ref. |T^. ( See in particular 
Fig.4 and Eq.(7) in Ref lH). For our Ptttt values and a ID expansion, = 1 independent 
of P^^. For the data in Table II (.r = 1.56. Thus V = (0.645 ± 0.130)dNjdr] fm^ and the 
range of V is 4.4±0.9 < V< 13.0±2.6 fm^ for 6.75 < dNjdr] < 20.2. 

(3) We assume that for dNc/drj > 6.75 the system is above the deconfinement transi- 
tion. The hot thermalized system expands, cools and then hadronizes. We attribute all of 
the measured volume to the expansion before hadronization. We neglect the subsequent 
expansion of the hadronic phase. Following Bjorken's derivation, we further assume that 
hydrodynamics of a massless relativistic ideal gas can describe the ID expansion and that 
the observed number of pions/fm^ are proportional to the entropy density s at hadronization. 
To estimate the pions/fm^ emitted by the source, the Bjorken ID boost invariant equation 
becomes 
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3/2 dNjdr] 
A 2 T 
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where A is the transverse area and T is the proper time |T^. The collisions occur at 
longitudinal coordinate z = and time t = 0. Eq. (2) describes an isentropic expansion 
s{r)/s{To) = %/r and 



(3) 



A 



where Tq is the initial proper time when thermahzation has occurred. For a relativistic 
massless ideal gas above the phase transition the maximum expansion velocity, responsible 
for most of the longitudinal expansion, is likely to be the sound velocity, = 1/3 The 



expansion time t = z/vg = in Rg/vs and T = (3z^ — z^Y^"^ = a/2 z. We note that 7} is the 
proper time at hadronization. 

Tf = V21rRg = V2 inh dNjdr] (4) 

and Eqn. (2) becomes 

3/2 dNjdri 1/V2 3/2 1/^2 



(5) 
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where 1 / \/2 is the effective Ar] slice. Thus is independent of dN^/ drj and one obtains 

= 1.64 ± 0.33 (stat) pions fm~^ (6) 

This value indicates that the deconfinement transition occurs at a definite entropy 
density. Since s oc is constant we can directly evaluate rin using the total number of pions 
emitted divided by the total volume for the data set in Table I. We choose the lowest P^^r 
value where r = 0.95 fm, Rq = 1.2 fm, and iji = 1.43. Here the average pseudo-rapidity 
density is < dNc/dr] >= 14.4 and Eqns. (2) and (4) become 

3/2 dNjdf] 1/V2 . 3 

= — - — = 1.57 ± 0.25 (stat) pions fm (7) 

TTT^ 2 tR Rg 

which has a smaller statistical error than (6). 

(4) The negative particle pt spectrum is used to measure the temperature. A slope 
parameter is obtained from a fit of the invariant cross section d'^Nc/dy d'^pt to the 
function Aexp{—bpt) for 0.15 < pt < 0.45 GeV/c |]T^. The value is constant to ± 1% for 
6.75 < dNc/df] < 20.2. Transverse flow has not been seen in p-p reactions at lower energies 



15| , p!6[| . In heavy ion reactions the transverse flow is attributed to final state interactions of 
the hadrons which presumably are not important in p-p collisions. The fact that Rq increases 
by a factor of three and remains constant to ± 1%, suggests that the transverse flow 
is negligible. The components a"2,o"3 in the NSD cross section indicate that the parton- 
parton mean free paths are shorter in high energy collisions. Since gluon-gluon interactions 
dominate in the initial encounters, early thermahzation ~ 0.5 fm/c when T ~ 200 MeV is 
likely §]. We interpret b~^ = T = 179.5 ± 5 MeV (syst) as the hadronization temperature. 
We neglect the expansion of the hadronic phase following hadronization i.e. decoupling 
is associated with hadronization. The systematic error estimate is based on possible kaon 
(K^) misidentification in the negative particle spectrum at low pt. We have not made a 
correction for the effect of resonance decays on the negative particle pt spectrum. We note 
that the negative particle temperature is significantly higher than the temperature based 
on the spectra of idenitified pions which include resonance particle decay pions (T~168 
MeV). A hadronization inverse slope parameter can be estimated from our measurement 
of the relative yields of mesons and hyperons as shown in Fig. 3, using all the events with 
dNc/drj > 6.75. The hadron yield versus rest mass inverse slope parameters indicates 162 < 



Tm < 173 MeV. Similar ~ 168 MeV values, based on thermal model analyses of liadron 
yield ratios, have been seen in high energy p p, p p, e^e~ and heavy ion reactions fl^. This 
has been interpreted as evidence for a universal limiting temperature Tm for hadrons, the 
Hagedorn temperature 
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FIG. 3. Relative meson and hyperon yields versus rest mass[19,20]. For the mesons, the inverse 
slope parameter = 162ib5 MeV, and for the hyperons Tm = 173ib 12 MeV. 

(5) We can also use the average measured energies to estimate the hadronization energy 



density ep |1T9|. Since ep = 3/4 T sp, ep is also constant versus dNc/dr] |T 



ep 



vrr^ X 2 X iji X h 



(8) 



where (mh) 



m'f^ + pIY^'^ is the average transverse mass of hadron h, Fh is a hadron 



abundance factor which also accounts for the neutral hadrons of each species. We have 
determined Fh for vr, i^, n. A", S, etc. For r = 0.95 fm, (.r = 1.56, h=0.073, ep 
becomes 



ep = 1.10 ± 0.26 (stat) GeV/fm^ 



(9) 



(6) We can estimate the average number ric of constituents in volume V at temperature 
T, for a system without boundaries [pl[] 



V 



G{T) 1.202 (kTf 



(10) 



where G{T) are the number of degrees of freedom (DOF). For a pion gas G{T) = 3, = 1 
fm^, and T = 179.5 MeV. The average number of pions (pion gas) in the source is n^r = 0.28 



pions/fm^. We observe 1.57 pions/fm^, emitted from the source at temperature T = 179.5 
MeV, which requires many more DOF. 

For a quark-gluon plasma G{T) = Gg{T) + Gq{T) + ^^(r) = 16 + 21/2 (/) where / are 
the number of quark flavors [|l^. Gg{T) are the gluon DOF; Gq(T), Gq(T) are the quark, 
antiquark DOF. We assume that pion emission from the source can be determined by the 
number of constituents in the source at hadronization, that one pion is a quark antiquark 
{q,q) pair and that two gluons {2g) are required to produce two pions (see insert Fig. 2). 

= Ug + {ug + nq)/2 (11) 



Our data indicates that ~ 6% strange quarks are present at hadronization p9|p0| . Thus 
we use / = 2 to evaluate Eqn. (10) where V = 1 fm^ and, 

n^ = {l + 2x 21/QA)Gg x 16.1 (GeV) (12) 

where Gg are the effective number of gluon DOF. For ra^ = 1.57/fm'^ and T = 0.1795 GeV, 
we obtain Gg = 10.18. The total number of DOF are, 

G{T) = ng + ng + Ug = (1 + 21/16)^3 = 23.5 ± 6 DOF (13) 

nearly eight times the DOF for a pion gas. A second method for estimating the DOF is to 
use the energy density and temperature at hadronization. For the isentropic expansion, the 



energy E in the volume V at temperature T is |21 



^^ymil^T'. (14) 

30 ^ ^ 

For ep = 1.10 ± 0.22 (stat) GeV/fm^ and T = 179.5 ± 5 MeV, we find G{T) = 24.8 ± 6.2 
(stat) quark gluon DOF, in good agreement with the DOF using the number of constituents 
(Eq. 13). 

(7) Two Lorentz contracted nucleons collide ai t = 0, z = and the thermalized con- 
stituents are assumed to emerge at Tq. Suppose we choose Tq = 1.0 fm/c. For a given expan- 
sion velocity, the data determines the hadronization proper time Tf and 1.09 < Tf < 3.25 
fm/c. For 6.75 < dNc/dr] < 20.2, the deconfined volumes V, determined by the data, range 
between AA < V < 13.0 fm^. For dNc/drj > 6.75 and using G{T) from Eq. (14) 

e/T^ = 7rV30 G{T) = 8.15 ± 2.0{stat) (15) 



in general agreement with lattice gauge calculations |^ . The ratio of the initial temperature 
Ti to the final Tf is T^/Tj = {Tf/ToY/^ and 185 < < 266 MeV. The ratio of the initial 
energy density to the final energy density ej is e^/e^ = {Tf /TqY/^ and is 1.23 < < 5.30 
GeV/fm'^ for 6.75 < dNc/dr] < 20.2. Note a different choice of % would change the Tj and 
ei estimates. 

In summary, the HBT analysis and the constant temperature versus dNc/dr] are consis- 
tant with a model in which a pion source undergoes a ID expansion with total longitudinal 
dimension 2 IrRg directly proportional to dNc/drj. We have used the Bjorken ID model to 
analyze our data. We find that there is a unique hadronization entropy density and temper- 
ature at which the pions are produced independent of dNc/drj. We have used phase space 



estimates of the average number of thermalized constituents in volume V at temperature T 
and the measured energy density ep to compute the number of DOF in the source. However, 
we note that reducing the average expansion velocity from = 1/3 to = 1/5 reduces 
the DOF estimate by 30%. Then the lower limit for the DOF is 16.6±4.2, still substantially 
larger than the pion gas DOF of 3. This lower hmit allows a more conservative argument 
that quark-gluon constituents are present in the large deconfined volumes. Our estimate 
of the number of DOF in the source (23.5±6, 24.8±6.2) is in general agreement with those 
expected for a quark-gluon plasma. The n^, ep, and T values characterize the quark-gluon 
to hadron thermal phase transition. We expect that these hadronization conditions will be 
observed in p-p coUisions at the CERN Large Hadron CoUider where higher pseudorapid- 
ity density dNc/dr) values will produce even larger deconfined volumes and longer plasma 
lifetimes. 

We would like to acknowledge the important support for Fermi National Laboratory 
experiment E-735 by J. D. Bjorken and L. M. Lederman. This work was supported in part 
by the United States Department of Energy and the National Science Foundation. 



REFERENCES 



* Current Address: Department of Physics, National Technical University of Athens, 
Athens, Greece. 

** Current Address: Department of Physics, University of Virginia, Charlottesville, Vir- 
ginia. 

[1] UAl collaboration, Phys. Lett. B 123 (1983) 115; Physics Letters 107B (1981) 320. 

[2] L. Van Hove, Phys. Lett. B 118 (1982) 138. 

[3] J. D. Bjorken, Fermilab Pub. 82/44-THY. 

[4] L. McLerran, Rev. Mod. Phys. 58 (1986) 1021. 

[5] E-735 proposal. Search for a deconfined quark-gluon plasma phase, Fermilab P-735 
(1983). 

[6] T. Alexopoulos, et al, (E-735 collaboration) Phys. Lett. B 435 (1998) 453; W. D. 

Walker, private communication. 
[7] X. N. Wang, Phys. Rep. 280 (1977) 287. 

[8] B. Miiller and A. Trayanov, Int. J. Mod. Phys. C5 (1994) 113; W. Poschel and B. 
Miiller, Phys. Rev. D 60 (1999) 114505. 

[9] T. Alexopoulos, et al, (E-735 collaboration) Phys. Rev. D 48 (1993) 1931. 
[10] E. V. Shuryak, Phys. Lett. 44B (1973) 387; G. Cocconi, Phys. Lett. 49B (1974) 459. 
[11] R. M. Weiner, Physics Reports 127 (2000) 249. 
[12] B. R. Schlei, et al, Phys. Lett. B 293 (1992) 275. 
[13] J. D. Bjorken, Phys. Rev. D 27 (1983) 140. 

[14] T. Alexopoulos, et al, (E-735 collaboration), Phys. Lett. B 336 (1994) 599. 

[15] A. T. Laasanen, et al., Phys. Rev. Lett. 38 (1977) 1. 

[16] I. Bearden, et al. (NA44 Collaboration). Phys. Rev. Lett. 78 (1977) 2080. 

[17] U. Heinz, Nucl. Phys. 685 (2001) 414c. 

[18] R. Hagedorn and J. Rafelski , Phys. Lett. 97B (1980) 136. 

[19] T. Alexopoulos, et al, (E-735 collaboration), Phys. Rev. D 46 (1992) 2773; Phys. Rev. 
D 48 (1993) 984. 

[20] T. Alexopoulos, et al, (E-735 collaboration), Z. Phys. C67 (1995) 411. 
[21] R. K. Pathria, Statistical Mechanics, Pergamon Press, Ltd, 1972, p 187f. 
[22] F. Karsch, et al, Phys. Lett. B 478 (2000) 249. 



TABLES 



TABLE I. Fitted values of radius Rq, lifetime r, and chaoticity A in the Gaussian parameter- 
ization with respect to qt and qq. Values are a function of average two-pion total momentum -Ptttt 
or average two-pion transverse momentum Pf. The total momentum interval containing the data 



is listed in column 1. Momentum is in GeV/c. The errors are statistical. 

P.n {P..) Rg (fm) r (fm) A (P,) 

0.2-0.5 0.404 1.20±0.05 0.95±0.06 0.24±0.01 0.369 

0.2-0.7 0.503 1.05±0.08 0.71±0.05 0.25±0.01 0.462 

0.5-1.0 0.708 0.80±0.07 0.67±0.07 0.23±0.02 0.650 

0.7-1.2 0.900 0.60±0.06 0.64±0.05 0.26±0.03 0.832 

0.9-1.7 1.175 0.58±0.06 0.53±0.07 0.26±0.02 1.087 

>1.0 1.403 0.48±0.06 0.45±0.05 0.21±0.02 1.285 

>1.2 1.600 0.43±0.06 0.41±0.06 0.23±0.02 1.479 



1 f\ 



TABLE II. Fitted values of the (longitudinal) radius Rq (transverse) lifetime r and corrected 
chaoticity A in the Gaussian parameterization with respect to qt and go- Values are a function of 
the average charged multiplicity per unit of pseudorapidity. Charged particle multiplicity intervals 
containing the data are listed in column 1. The errors are statistical. 



A^c {dNjdr]) RGifm) r (fm) A 

0-60 6.75 0.62 ± 0.09 0.53 ± 0.07 0.39 ± 0.05 

0-80 9.00 0.70 ± 0.09 0.65 ± 0.06 0.32 ± 0.03 

60-100 12.5 1.00 ± 0.08 0.86 ± 0.12 0.25 ± 0.01 

80-120 15.5 1.10 ± 0.10 0.89 ± 0.12 0.23 ± 0.01 

100-240 18.2 1.52 ± 0.14 0.99 ± 0.15 0.21 ± 0.02 

120-240 20.17 1.86 ± 0.35 0.88 ± 0.20 0.19 ± 0.03 
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